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We report on the first step of a systematic study of how gaugino mass unification can be probed at the LHC 
in a quasi-model independent manner. Here we focus our attention on the theoretically well-motivated mirage 
pattern of gaugino masses, a one-parameter family of models of which universal (high scale) gaugino masses are 
a limiting case. Using a statistical method to optimize our signature selection we arrive at three ensembles of 
observables targeted at the physics of the gaugino sector, allowing for a determination of this non-universality 
parameter without reconstructing individual mass eigenvalues or the soft supersymmetry-breaking gaugino masses 
themselves. In this controlled environment we find that approximately 80% of the supersymmetric parameter space 
would give rise to a model for which our method will detect non-universality in the gaugino mass sector at the 
10% level with 0(10 fb -1 ) of integrated luminosity. 



1. MOTIVATION 

Given that the LHC era is nearly here, it is 
hardly surprising to hnd members of the high 
energy theory community are increasingly turn- 
ing their attention away from esoteric issues of 
model-building and towards issues of directly con- 
fronting theory with data. This activity has pro- 
duced a number of new measurement techniques. 
Previously unconsidered signatures have been de- 
veloped by studying certain "what-if" scenarios, 
and a new emphasis has been placed on mak- 
ing concrete predictions from famously nebulous 
string constructions. All of this represents impor- 
tant progress, but here I wish to consider another 
aspect of the challenge that lies ahead. Specif- 
ically, I wish to imagine the state of our held 
three to five years from now. Assuming some- 
thing like the minimal supersymmetric standard 
model (MSSM) is discovered, as so many expect, 
then we can imagine being presented with a num- 
ber of measurements that are indirectly related 
to the underlying soft supersymmetry-breaking 
mechanisms. Here I will discuss what is likely 
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just the first step in a research program into 
what comes next: how to connect the multi- 
ple LHC observations to organizing principles in 
some (high-energy) effective Lagrangian of under- 
lying physics. 

In addressing this issue we must be careful to 
avoid the pitfalls of the "LHC Inverse" problem: 
the empirical fact that even in very restrictive 
model frameworks it is quite likely that more 
than one set of model parameters will give pre- 
dictions for LHC observations that are in good 
agreement with the experimental data [2], Much 
recent work has focused on how to address this 
issue |3l4l5l6l7t8j . and we will borrow much of 
the philosophy and many of the useful techniques 
from this recent literature. The LHC inver- 
sion problem generally results from trying to use 
a large ensemble of (correlated) observations to 
constrain a multi-dimensional parameter space. 
But our ultimate goal as theorists is to under- 
stand broad properties of the underlying physics 
itself. Ever more precise measurement of a single 
parameter - say the gluino mass - does not neces- 
sarily further this goal. The most important such 
broad property is undoubtedly the issue of gau- 
gino mass unification: few properties of the super- 
partner spectrum have more far-reaching implica- 
tions for low-energy phenomenology, the nature of 
supersymmetry breaking, and the structure of the 
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underlying physics Lagrangian |9ll0j . But these 
soft parameters are not themselves directly mea- 
surable at the LHC [TT] . One might consider per- 
forming a fit to some particular theory, such as 
minimal supergravity (mSUGRA), in which uni- 
versal gaugino masses are assumed [12] . but we 
are not so much interested in whether mSUGRA 
- or any other particular theory for which gaugino 
mass universality is a feature - is a good fit to the 
data. Rather, we wish to know whether gaugino 
mass universality is a property of the underlying 
physics independent of all other properties of the 
model. 

We will therefore begin our attack by con- 
sidering a concrete parametrization of non- 
universalities in soft gaugino masses. In recent 
work by Choi and Nilles |10j soft supersymmetry- 
breaking gaugino mass patterns were explored in 
a variety of string-motivated contexts. In par- 
ticular, the so-called "mirage pattern" of gau- 
gino masses [13114115] provides an interesting case 
study in gaugino mass non-universality. This 
paradigm gets its name from the fact that should 
the mirage pattern of gaugino masses be used as 
the low-energy boundary condition of the (one- 
loop) renormalization group equations then there 
will exist some high energy scale at which all 
three gaugino masses are identical [16] . The set of 
all such low-energy boundary conditions that sat- 
isfy the mirage condition defines a one-parameter 
family of models. In the parametrization we 
adopt from [ID] , the gaugino mass ratios at the 
electroweak scale take the form 

Mi : M 2 : M 3 ~ (l+0.66a) : (2+0.2a) : (6-1.8a)(l) 

where the case a = is precisely the unified 
mSUGRA limit. In the limit of very large values 
for the parameter a the ratios among the gaugino 
masses approach those of the anomaly-mediated 
supersymmetry breaking (AMSB) paradigm [17[ 
118] . In fact, the mirage pattern is most naturally 
realized in scenarios in which a common contribu- 
tion to all gaugino masses is balanced against an 
equally sizable contribution proportional to the 
beta-function coefficients of the three Standard 
Model gauge groups. Such an outcome arises 
in string-motivated contexts, such as KKLT-type 
moduli stabilization in D-brane models [19120] 



and Kahler stabilization in heterotic string mod- 
els [21|22|23|24|25|26l27|28j . Importantly, how- 
ever, it can arise in non-stringy models, such as 
deflected anomaly mediation [29 30] . We note 
that in none of these cases is the pure- AMSB limit 
likely to be obtained, so our focus here will be 
on small to moderate values of the parameter a 
in (P). 

2. METHOD 

The ultimate goal is to ask whether or not soft 
supersymmetry breaking gaugino masses obey 
some sort of universality condition independent of 
all other facts about the supersymmetric model. 
We begin by asking a simpler question: assum- 
ing the world is defined by the MSSM with gau- 
gino masses obeying the relation fl} , how well can 
we determine the value of the parameter a. At 
the very least we would like to be able to estab- 
lish that a with a relatively small amount 
of integrated luminosity. The first step in such 
an incremental approach is to demonstrate that 
some set of "targeted observables" [31] (we will 
call them "signatures" in what follows) is sensi- 
tive to small changes in the value of the parame- 
ter a in a world where all other parameters which 
define the SUSY model are kept fixed. 

2.1. Setting Up the Problem 

We will construct what we will call a "model 
line" by specifying the supersymmetric model in 
all aspects other than the gaugino sector. We 
choose a simplified set of 17 input parameters 
given by 

tan/3, m 2 H , m 2 H "1 
M 3 , A t , A b , Ar I (2) 

TO Qi,2> m U!,2, m Dl 2 , m Ll 2 , m El 2 (' 
m Q 3 j m U 3 , m D3 , m L3 , m E3 J 

which are understood to be taken at the elec- 
troweak scale (specifically A EW = 1000 GeV) so 
no renormalization group evolution is required. 
The gluino soft mass M3 will set the overall scale 
for the gaugino mass sector. The other two gau- 
gino masses M\ and M2 are then determined rel- 
ative to M3 via fl]). A model line will take the 
inputs of ^ and then construct a family of theo- 
ries by varying the parameter a from a = (the 
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mSUGRA limit) to some non-zero value in even 
increments. 

For each point along the model line we passed 
the model parameters to PYTHIA 6.4 for spec- 
trum calculation and event generation. Events 
are then sent to the PGS4 package to simulate 
the detector response. Additional details of the 
analysis will be presented in later sections. The 
end result of our procedure is a set of observable 
quantities that have been designed and (at least 
crudely) optimized so as to be effective at sepa- 
rating a = from other points along the model 
line in the least amount of integrated luminosity 
possible. 

2.2. Distinguishability 

The technique we employ to distinguish be- 
tween candidate theories using LHC observables 
involves the construction of a variable similar to 
a traditional chi-square statistic 
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We can identify any signature Ni with an ef- 
fective cross section tr, via the relation 



Ni/L, 



(4) 



where L is the integrated luminosity. We refer to 
this as an effective cross-section as it is defined 
by the counting signature Ni which contains in 
its definition such things as the geometric cuts 
that are performed on the data, the detector ef- 
ficiencies, and so forth. Furthermore these effec- 
tive cross sections, whether inferred from actual 
data or simulated data, are subject to statisti- 
cal fluctuations. The transformation in ^ allows 
for a comparison of two signatures with differing 
amounts of integrated luminosity. This will prove 
useful in cases where the experimental data is pre- 
sented after a limited amount of integrated lumi- 
nosity La, but the simulation being compared to 
the data involves a much higher integrated lumi- 
nosity Lb ■ We will assume that the errors associ- 
ated with the signatures Ni are purely statistical 
in nature and that the integrated luminosities La 
and Lb are precisely known, so that (ASab) 2 is 
given by 



where S is some observable quantity (or signa- 
ture). The index i — l,...,n labels these sig- 
natures, with n being the total number of signa- 
tures considered. The labels A and B indicate 
two distinct theories which give rise to the signa- 
ture sets Sf and Sf, respectively. Finally, the 
error term SS^ B is an appropriately-constructed 
measure of the uncertainty of the term in the nu- 
merator, i.e. the difference between the signa- 
tures. In this work we will always define a signa- 
ture S as an observation interpreted as a count 
(or number) and denote it with capital N. One 
example is the number of same-sign, same-flavor 
lepton pairs in a certain amount of integrated lu- 
minosity. Another example is taking the invariant 
mass of all such pairs and forming a histogram of 
the results, then integrating from some minimum 
value to some maximum value to obtain a num- 
ber. In principle there can be an infinite number 
of signatures defined in this manner. In practice 
experimentalists will consider a finite number and 
many such signatures are redundant. 
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where each cross section includes the (common) 
Standard Model background, i.e. <7j = 5f USY + 

CT SM . 

The variable (ASab) 2 forms a measure of the 
distance between any two theories in the space 
of signatures defined by the Si- We can use this 
metric on signature space to answer the following 
question: how far apart should two sets of sig- 
natures Sf and SP be before we conclude that 
theories A and B are truly distinct? To answer 
this question we note that in the limit in which 
the luminosities La and Lb are large the proba- 
bility distribution for the quantity (ASab) 2 given 
by " 



P(AS 2 ) = nx 2 n , x (nAS 2 ), 



(6) 



where Xn \ is the non-central chi-squared dis- 
tribution for n degrees of freedom. The non- 
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centrality parameter A is given by 



A = 



fa 



r B\2 



af/LA + af/Lt 



(7) 



and now the Cj represent ezact cross sections. 
From (O and (J7J) it is apparent that all the physics 
behind the distribution of possible (ASab) 2 val- 
ues is contained in the values of A and n. 

The problem of defining "distinguishability" 
has now been reduced to a well-understood prob- 
lem in statistics which can be solved analyti- 
cally. To say that two potential models have 
been distinguished - or that a set of n obser- 
vations are inconsistent with those derived from 
a simulation - we first demand that the quan- 
tity (ASab) 2 constructed from the n observa- 
tions satisfy (ASab) 2 > (ASAA) 2 \ 95th - In other 
words, we first require that the distance in signa- 
ture space be larger than that expected simply by 
quantum fluctuations with 95% confidence. The 
values of (AS^) 2 | g5th for A = are easily com- 
puted and tabulated for any value of n. 

2.3. Optimization 

Now return to the basic problem: using experi- 
mental data to distinguish two models that truly 
are distinct. Though we have A ^ 0, there is 
always a finite chance that an experimental mea- 
surement will not reveal this fact due to quan- 
tum fluctuations. For any given value of A ^ 0, 
the probability that a measurement of (ASab) 2 
will fluctuate to a value so small that it is not 
possible to separate two distinct models (to con- 
fidence level p) is simply the fraction of the prob- 
ability distribution in ((6| that lies to the left of 
the value (ASAA) 2 \ g5th - There is always some 
minimum value of the non-centrality parameter 
that can be chosen so that this fraction is be- 
low some pre-determined value. Let us call that 
value X min (n,p). Again, these values are simply 
determined by integrating the probability distri- 
butions (O and are independent of the physics of 
the problem. So, for example, given two distinct 
models A and B, any combination of n experi- 
mental signatures such that A > A tI1 j n (n,p = 0.95) 
will be effective in demonstrating that the two 
models are indeed different 95% of the time, with 
a confidence level of 95%. 



Let us assume for the moment that "model A" 
is the experimental data, which corresponds to 
an integrated luminosity of L cxp . Our "model B" 
can then be a simulation with integrated lumi- 
nosity L slm — qL cxp . Define the quantity 
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where Rab has the units of a cross section. Our 
condition for 95% certainty that we will be able 
to separate two truly distinct models at the 95% 
confidence level becomes 
A min (n, 0.95) 



Rab 



(9) 



Given two models A and B and a selection of 
n signatures both A m i n (n, 0.95) and Rab are 
completely determined. Therefore the minimum 
amount of integrated luminosity needed to sepa- 
rate the models experimentally will be given by 

A min (n,_p) 



Rai 



(10) 



Equation (|10[) provides us with a quantitative 
measure of the efficacy of any particular set of n 
(uncorrelated) signatures. We therefore wish to 
select a set of n signatures Si such that the quan- 
tity £ m in(p) as defined in (fT0|) . for a given value of 
p, is as small as it can possibly be over the widest 
possible array of model pairs A and B. We must 
also do our best to ensure that the n signatures 
we choose to consider are reasonably uncorrelated 
with one another so that the statistical treatment 
of the preceding section is applicable. 

The optimal strategy is generally to choose a 
rather small subset of the total signatures one 
could imagine. In part this is because in any 
set of observations the probability that any two 
are highly correlated with one another will only 
grow with the size of the set. But in addition 
not all possible signatures are equally effective at 
separating models. An absolute measure of the 
"power" of any given signature to separate two 
models A and B can be provided by considering 
the condition in (TlT)|) . For any signature Si we 
can define an individual (L m - nl )i by 



(-^min)i — A m ; n (l,p) 



9 l _ 



T B\2 



(11) 
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where, for example, A min (l,0.95) = 12.99. This 
quantity is exactly the integrated luminosity re- 
quired to separate models A and B, to confidence 
level p, by using the single observable Si . As more 
signatures are added, the threshold for adding the 
next signature in the list gets steadily stronger. 
Mathematically, the ratio A m i n (n,p)/A m i n (l, p) 
grows with increasing n. This indicates that as 
we add signatures with ever diminishing (£ m i n )i 
values we will eventually encounter a point of neg- 
ative returns, where the resulting overall L m i n 
starts to grow again. 

For a particular pair of models, A and B, it is 
always possible to find the optimal list of signa- 
tures from among a given grand set by ordering 
the resulting (£ m i n ), values and adding them se- 
quentially until a minimum of L mm is observed. 
To do so, we note that kinematic distributions 
must be converted into counts (and all counts 
are then converted into effective cross sections). 
This conversion requires specifying an integration 
range for each histogram. The choice of this range 
can itself be optimized, by considering each inte- 
gration range as a separate signature and choos- 
ing the values such that (L m i n )i is minimized. By 
repeating this optimization procedure over many 
background model lines it is possible to construct 
a list that will be at least close to optimal over a 
wide range of supersymmetric model space. 

3. SIMULATION AND SIGNATURE SE- 
LECTION 

To perform the optimization procedure we con- 
structed a large number of model families in the 
manner described in Section 12. 1( each involving 
the range —0.5 < a < 1.0 for the parameter a 
in steps of Aa = 0.05. For each point along 
these model lines we generated 100,000 events us- 
ing PYTHIA 6 . 4 and PGS4 using the default level- 
one triggers. To this we added an appropriately- 
weighted Standard Model background sample 
consisting of 5 fb _1 each of t/i and b/b pair pro- 
duction, high-pT QCD dijet production, single 
and Z-boson production, pair production of 
electroweak gauge bosons (W + W~, W ± Z and 
Z Z\ and Drell-Yan processes. Further object- 
level cuts were then performed, followed by an 



event-level cut on the surviving detector objects 
similar to those used in [2]. Specifically we re- 
quired all events to have missing transverse en- 
ergy $t > 150 GeV, transverse sphericity St > 
0.1, and H T > 600 GeV (400 GeV for events with 
2 or more leptons) where Ht =$t + X) lots Pr 1 

To examine the effectiveness of our candidate 
signature sets at measuring the value of the pa- 
rameter a we fixed "model A" to be the point 
on each of the model lines with a = and then 
treated each point along the line with a ^ as a 
candidate "model B." Clearly each model line we 
investigated - and each a value along that line - 
gave slightly different sets of maximally effective 
signatures. The lists we will present represent an 
ensemble average over these model lines. 

As a straw man we may consider the sin- 
gle most effective signature at separating models 
with different values of the parameter a. It is the 
effective mass formed from all objects in the event 

M™y=$ T + J2pf, as) 

all 

where we form the distribution from all events 
which pass our initial cuts. To turn this distri- 
bution into a count we integrate from M^ y — 
1250 GeV to the end of the distribution. We will 
refer to the variable in (fT2]) as signature "list" A. 
That this one signature would be the most pow- 
erful is not a surprise given the way we have set 
up the problem. It is the most inclusive possible 
signature one can imagine (apart from the over- 
all event rate itself) and therefore has the largest 
overall cross-section. Furthermore, the variable 
in (fT2|) is sensitive to the mass differences between 
the gluino and the lighter electroweak gauginos 
- precisely the quantity that is governed by the 
parameter a. Yet as we will see in Section [4] this 
one signature can often fail to be effective at all in 
certain circumstances, resulting in a rather large 
required L m { n to be able to separate a = from 
non- vanishing cases. In addition it is built from 
precisely the detector objects that suffer the most 
from experimental uncertainty. This suggests a 
larger and more varied set of signatures would be 
preferable. 

We next consider the five signatures in Table [TJ 
These signatures were chosen by taking our most 
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Table 1 



Signature List B. Distributions are integrated from "Min Value" to "Max Value". 





Description 


Min Value 


Max Value 


1 


M^ s [0 leptons, > 5 jets] 


1100 GeV 


End 


2 


M™ y [0 leptons, < 4 jets] 


1450 GeV 


End 


3 


M^ y [> 1 leptons, < 4 jets] 


1550 GeV 


End 


4 


pr (Hardest Lepton) [> 1 lepton, > 5 jets] 


150 GeV 


End 


5 


M^l s [0 leptons, < 4 jets] 


GeV 


850 GeV 



effective observables and restricting ourselves to 
that set for which e = 10%. We again see the 
totally inclusive effective mass variable of (fT2)) as 
well as the more traditional effective mass vari- 
able, M 3 c g S , defined via (fl~2|) but with the scalar 
sum of pt values now running over the jets only. 
We now include the px of the hardest lepton in 
events with at least one lepton and five or more 
jets, as well as the invariant mass of the jets 

in events with zero leptons and 4 or less jets. The 
various jet-based effective mass variables would 
normally be highly correlated with one another if 
we were not forming them from disjoint partitions 
of the overall data set. The favoring of jet-based 
observables to those based on leptons is again 
largely due to the fact that jet-based signatures 
will have larger effective cross-sections for reason- 
able values of the SUSY parameters in ^ than 
leptonic signatures. The best signatures are those 
which track the narrowing gap between the gluino 
mass and the electroweak gauginos and the nar- 
rowing gap between the lightest chargino/second- 
lightest neutralino mass and the LSP mass. In 
this case the first leptonic signature to appear - 
the transverse momentum of the leading lepton 
in events with at least one lepton - is an example 
of just such a signature. 

Finally, let us consider the larger ensemble of 
signatures in Table [H In this final set we have 
relaxed our concern over the issue of correlated 
signatures, allowing as much as 30% correlation 
between any two signatures in the list. This al- 
lows for a larger number as well as a wider variety 
of observables to be included. This can be very 
important in some cases in which the supersym- 
mctric model has unusual properties, or in cases 



where the two a values being considered give rise 
to different mass orderings (or hierarchies) in the 
superpartner spectrum. In displaying the signa- 
tures in Table [2] we find it convenient to group 
them according to the partition of the data being 
considered. Note that the counting signatures are 
taken over the entire data set. 

The first counting signature is simply the to- 
tal size of the partition in which the events have 
at least one lepton and 4 or less jets. The next 
two signatures are related to "spoiler" modes for 
the trilepton signal. Note that the trilepton sig- 
nal itself did not make the list: this is a wonderful 
discovery mode for supersymmetry, but the event 
rates between a model with a = and one with 
non- vanishing a were always very similar (and 
low) . This made the trilepton counting signature 
ineffective at distinguishing between models. By 
contrast, counting the number of b-jet pairs (a 
proxy for counting on-shell Higgs bosons) or the 
number of opposite-sign electron or muon pairs 
whose invariant mass was within 5 GeV of the 
Z-mass (a proxy for counting on-shell Z-bosons) 
were excellent signatures for separating models 
from time to time. This was especially true when 
the two models in question had very different val- 
ues of a such that the mass differences between 
N2 and Ni were quite different in the two cases. 
Signature #8 is defined as the following ratio 



'jet 



ict3 1 iet4 

Pt + F T 

ictl . ict2 

Ft + Pt 



(13) 



where pl^ 1 is the transverse momentum of the i- 
th hardest jet in the event. For this signature 
we require that there be at least three jets with 
Pt > 200 GeV. This signature, like the pt of the 
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Table 2 



Signature List C. Distributions are integrated from "Min Value" to "Max Value". 





Description 


Min Value 


Max Value 


Counting Signatures 


1 


N t [> 1 leptons, < 4 jets] 




2 


N i+e - [Af£T = M z ± 5 GeV] 




3 


N B [> 2 B-jets] 






[0 leptons, < 4 jets] 


4 


M ™y 


1000 GeV 


End 


5 




750 GeV 


End 


6 




500 GeV 


End 


[0 leptons, > 5 jets] 


7 


litany 


1250 GeV 


3500 GeV 


8 


r je t [3 jets > 200 GeV] 


0.25 


1.0 


9 


PT(4th Hardest Jet) 


125 GeV 


End 


10 




0.0 


0.25 


[> 1 leptons, > 5 jets] 


11 




0.0 


0.25 


12 


(Hardest Lepton) 


150 GeV 


End 


13 


pr(4th Hardest Jet) 


125 GeV 


End 


14 


# T + M^ S 


1250 GeV 


End 



hardest lepton or the pr of the 4th hardest jet, 
was effective at capturing the increasing softness 
of the products of cascade decays as the value of 
a was increased away from a = 0. 

4. RESULTS 

In this section we examine how well our sig- 
nature lists in perform in measuring the value of 
the parameter a which appears in (|TJ) . These lists 
were designed through the analysis of several hun- 
dred "model lines" as described previously. To 
test the approach we generated a random collec- 
tion of 500 new models with six points on the 
a-lines ranging from to 0.5. In this case a 16- 
dimensional parameter space defined by the quan- 
tities in {2]| was considered. Specifically, slepton 
and squark masses were allowed to vary in the 
range 300 GeV to 1200 GeV with the masses of 
the first and second generation scalars kept equal. 
The gaugino mass scale given by M3 and the jj,- 



parameter were also allowed to vary in this range. 
The pseudoscalar Higgs mass tua was fixed to be 
850 GeV and the value of tan/3 was allowed to 
vary from 2 to 50. If all points along the a-line 
satisfied all experimental constraints on the su- 
perpartner mass spectrum, then 100,000 events 
were generated for each of the six points along 
the a-line in the manner described in Section [5J 
Using this data the value of L m - m was computing 
using (fTU|) for each of our three signature sets. 

The results of this analysis are given in Fig- 
ures Q] and [2] In Figure Q] the integrated luminos- 
ity needed to detect a ^ for 95% of our ran- 
dom models is given as a function of the five non- 
vanishing a values simulated. Since the random 
model sample includes examples with very differ- 
ent superpartner mass scales, the overall super- 
symmetric production cross-section varies much 
more across this sample than in the controlled 
model sample described above. We therefore take 
this into account by plotting the required number 



1 

0.0 0.1 0.2 0.3 0.4 0.5 

Gaugino Mass Non-Universality Parameter a 



Fi gure 1. -Lmin required to detect a ^ for 95% of 
the random models. 



of supersymmetric events in Figure [5] While our 
single best signature ("List A") has some resolv- 
ing power, it is the larger sets that will prove most 
effective in establish the non-universality of gau- 
gino masses, as to be expected. Yet going from 
the 5 signatures of List B to the 14 signatures of 
List C produces only a marginal increase in re- 
solving power. This is evidence of the increased 
correlations amongst the signatures of List C. It 
is also demonstrative of the general principle of 
diminishing returns to adding new observables to 
any particular set. Even using our best set of sig- 
natures (List C) it will require nearly 100 fb _1 to 
be able to detect non-universality at the level of 
a ~ 0.1 for an arbitrary supersymmetric model. 
Yet for the vast majority of models the departure 
from universality should become apparent after 
just 10-20 fb" 1 . Departures from universality at 
the level of a ~ 0.3 should be apparent using 
this method for most supersymmetric models af- 
ter just a few fb _1 . 

To understand why this approach works, it 
is useful to examine the signature results them- 
selves. As an example, consider a pair of signa- 
tures drawn from List C in Table [U Figure [3] 
shows a two-dimensional slice of the signature 
space "footprint" for our large set of model vari- 
ations |32l33l34j . In these figures the results have 
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Figure 2. iV m i n required to detect a ^ for 95% of 
the random models. 



been normalized to 5 fb _1 of data. The count rate 
for signature #9 is shown versus that of signature 
#11 of List C for the case a = and a — 0.33 
(top panel), a — 0.66 (middle panel) and a = 1 
(bottom panel) . We have chosen this pair for the 
dramatic separation that can be achieved, though 
similar results can be obtained with other pairs 
of signatures. 

The power of our inclusive signature list ap- 
proach lies in the choice of signatures and their 
ability to remain highly sensitive to changes in the 
physical behavior of each model. This feature is 
reflected qualitatively in the visual clustering of 
the data points, which become progressively more 
distinct as the parameter a is increased. As the 
regions separate it becomes increasingly less likely 
that a model from one class can be confused with 
a model from the other class, even when consider- 
ing statistical fluctuations. In our approach this 
manifests itself when one computes Rab, which 
reflects the "distance" in signature space between 
the two models under comparison, and which be- 
comes large when the models are sufficiently dif- 
ferent from one another. 

5. CONCLUSIONS 

If supersymmetry is discovered at the LHC the 
high energy community will be blessed with a 
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Figure 3. Footprint-style plot for a pair of signa- 
tures from List C. Total counts for signature #11 versus 
signature #13 of List C is given for the case a = (green 
triangles) a # (black squares). The cases shown are 
for a = versus a = 0.33 (top panel), a = 0.66 (middle 
panel) and a = 1 (bottom panel). The axes measure the 
number of events for which the kinematic quantity was 
in the range given in Tabic [2] Larger values of the non- 
universality parameter a correspond to a greater degree 
of separation between the two model "footprints." 



large number of new superpartners whose masses 
and interactions will need to be measured. Un- 
doubtedly performing global fits of the many ob- 
servables to the parameter space of certain priv- 
ileged and well-defined benchmark models will 
be of great help in making sense of this embar- 
rassment of richness. It is an interesting ques- 
tion to ask whether it is possible to fit to cer- 
tain model characteristics rather than to any par- 
ticular model itself. Perhaps the most impor- 
tant such characteristic is the pattern of soft 
supersymmetry-breaking gaugino masses. No 
other property of the low-energy soft Lagrangian 
is more easily linked to underlying high-scale 
physics, particularly if that high-scale physics is 
of a string-theoretic origin. We are thus inter- 
ested in asking whether we can identify the pres- 
ence on non-universalities in the gaugino sector 
independent of all other properties of the su- 
perpartner spectrum. In the present work we 
have decided to begin this process with a sim- 
ple parametrization in terms of model "lines" , in 
the spirit of previous benchmark studies such as 
the Snowmass Points & Slopes [35] in which only 
the single non-universality parameter is varied. 
By understanding how the observable physics at 
the LHC is affected by this parameter - and then 
repeating the analysis many times with the other 
supersymmetric parameters varied - we were able 
to obtain two sets of observables that fared well in 
detecting the presence of non-universalities with 
relatively small amounts of integrated luminosity. 
The larger of the two sets generally performed 
slightly better, but at the expense of allowing 
signatures that are correlated with one another 
at the 30% level. Requiring a correlation at only 
the 10% level (and thus using a shorter list of 
observables) had only a small effect on the abil- 
ity to distinguish the size of the non-universality 
parameter a. Broadly speaking, we find that a 
non-universality at the 10% level can be measured 
with 10-20 fb _1 of integrated luminosity over ap- 
proximately 80% of the supersymmetric parame- 
ter space relevant for LHC observables. If we are 
interested in measurements at only the 30% level 
these numbers change to 5-10 fb^ 1 over approxi- 
mately 95% of the relevant parameter space. 
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